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ABSTRACT 

In this paper we investigate the power spectrum of the unresolved 0.5-2 keV CXB with 
deep Chandra 4 Ms observations in the CDFS. We measured a signal which, on scales 
>30", is significantly higher than the Shot- Noise and is increasing with the angu- 
lar scale. We interpreted this signal as the joint contribution of clustered undetected 
sources like AGN, Galaxies and Inter-Galactic- Medium (IGM). The power of unre- 
solved cosmic sources fluctuations accounts for ~12% of the 0.5-2 keV extragalactic 
CXB. Overall, our modeling predicts that ^20% of the unresolved CXB flux is made 
by low luminosity AGN, ^25% by galaxies and ~55% by the IGM (Inter Galactic 
Medium). We do not find any direct evidence of the so called Warm Hot Intergalac- 
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tic Medium (i.e. matter with 10 5 K<T<10 7 K and density contrast 5 <1000), but we 
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estimated that it could produce about 1/7 of the unresolved CXB. 

We placed an upper limit to the space density of postulated X-ray-emitting early 
black hole at z>7.5 and compared it with SMBH evolution models. 

Key words: (cosmology:) dark matter, (cosmology:) large-scale structure of universe, 
X-rays: galaxies, galaxies: active, (cosmology:) diffuse radiation 



1 INTRODUCTION 

The Cosmic X-ray Background (CXB) is the result of a 
multitude of energetic phenomena occurring in the Universe 
since the epoch of the formation of the first galaxies. Its na- 
ture has been investigated in the last 50 years with several 
telescopes but only in t he 80's it became clear that its main 
contributors are AGN (|Setti fc Woltierlll989l ). Later it has 
been found that also galaxies, galaxy clusters, large scale 
structures and diffuse hot gas in the Milky Way ar e sources 
contributing to the CXB (|Fabian fc Barcons|[l99^ . 
With the launch of ROSAT, Chandra and XMM- Newton, 
the knowledge about the nature of sources contributing to 
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the flux of the CXB became suddenly clear . Deep surveys 
like the Chandra Deep Field South (C DFS, iGiacconi et al.l 
| 200ll:lLuo et al.ll2008l : IXue et alj|20ljh and North (CDFN 
[ Brandt et al. 200 ll ). and the Lockman Hole (|Brunner et al.l 
have almost conclusively resolved the problem of the 
CXB below 10 keV. In fact at the fl ux limits of Chan- 
dra and XMM-Newton, about 90-95% ifBauer et alj 120041 : 
iMoretti et ail 120031 : iLehmer et ail |2012j ) of the 0.5-2 keV 
CXB flux has been resolvecQ into point and extended 

sources. 

CXB synthesis models (see e.g. iTreister fc Urrvl 120061 : 



1 In these papers, the fraction of unresolved CXB has been esti- 
mated at the flux of the faintest detected source, here we measure 
the average value on the investigated area 
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iGilli et al.|[2007h predict that at fluxes fainter that the cur- 
rent limits, the AGN source counts progressively flatten and 
then galaxies become the more abundant sources. 
Galaxies emerge as the dominant population at faint 0.5- 
2 keV fluxes, with b reak-even point at around ~10~ 17 
erg cm" 2 s'^see a lso iLehmer et ai1l2007l ; IXue et al.ll201ll ; 
iLehmer et ai]|2012r i and can make up for a sizable fraction 
of the X-ray background (9-16%; Ranalli et al. 2005). The 
contribution from galaxy clusters, down to a mass limit of 
~10 13 Mq, has been estimated to be of the order of ~10% 
of th e total 0.5-2 keV CXB (|Gilli et al.|[l999l ; iLemze et all 
l2009h . 

Several authors studied the clustering properties of the 
CXB to unveil the nature of the sources produc- 
ing its flux and their prope rt ies jM artin-Mirone s et al] 
19911: IWu fc. Anderson! 1 19921 : ISoltan fe Hasingerl 1 1991 : 
Barcons et al. 19981 ). The most used technique is that of the 



two-point autocorrelation function of the surface brightness 
of the CXB. While at the time of HEAO-1 the task was to 
determine what was the contribution of QSO to the CXB, af- 
ter ROSAT most of the investigations in this field have been 
performed t o unveil signatures of the Warm Hot Intergalac- 
tic M edium (jSoltan et al]|2002l : ISoltanll2006l; iGaleazzi et"ahl 
I2009T ) and to test Cosmological models I Diego et al.ll2003h . 
This kind of studies is therefore extremely powerful to study 
population of sources which are beyond the resolving and 
detection capabilities of instruments. 

In the GOODS fields, iHickox fc Markevitchl l|2006l . 
12007ft have shown that, after excising detected point and 
extended sources plus faint HST detected galaxies, the 
spectrum of the soft CXB was still showing a signal in 
the 0.5-2 keV band. While a large amount of it could 
be attributed to a local component (Solar Wind Charge 
Exchange and Milky Way thermal emission), they have 
shown that below 1 keV the fraction of the resolved CXB is 
not sensitive to the removal of HST galaxies, thus arguing 
for a purely diffuse nature of the unresolved CXB. In the 
1-2 keV band they have shown that, by removing HST 
sources areas from the analysis, the fraction of unresolved 
CXB drops by ~50%. This means that a large fraction of 
the unresolved 1-2 keV CXB flux could be produced in 
faint undetected galaxies (active or non-active). Moreover, 
they speculated that a large fraction of the unresolved 
CXB may be due to faint 3.6/im IRAC sources, suggesting 
that high-z or absorbed sources could produce a large 
part of such a radiation. According to their estimate, the 
remaining fraction of the CXB (i.e. 2-3% of the 0.5-2 
keV CXB) remains consistent with the prediction of 
Warm Hot Intracluster Medium ( WHIM) intensity from 
hydrodynamical simulations (see e.g Cen fc Ostrikerl 1 19991 : 
lUrsino fc Galeazzill2006l : iRoncarelli et al.ll2006l . l2012l ). 
In the local Universe about 30-40% of the b aryons are miss- 
ing w ith respect to what is measured at z~3 (|Fukugita et al] 
1 19981 ). Simulations predict that most of these baryons got 
shock- heated and at z<l they should have a temperature 
of the order 10 5 - 10 7 K and therefore emit thermal X-rays 
l|Cen fc Ostrikerl Il999t ). Controversial evidences on the 
prop erties of such a medium have been published so far (see 
e.g., iNicastro et al.ll2005l : iKaastra et all l200d : iFang et al] 
l20ld : IShull et al.l I2011T K Though with very low surface 
brightness, the WHIM can be distinguished from other 
kinds of diffuse emissions on the basis of its clustering 



properties (|Ursino fc Ga leazzi 2006). In fact, its emission 
follows that of clusters and filaments and peaks at low red- 
shift (z~0.5), thus showing a typical feature in the angular 
clustering of unresolved CXB fluctuations. However, WHIM 
is not the only expected component of the unresolved CXB 
arising from thermal emission of the Inter Galactic Medium 
(IGM). X-ray surveys in the local Universe revealed X-ray 
emission from local galax y groups down to ma sses of the 
order of 10 12 M (see e.g. lEckmiller et al.llioill 'l. Since the 
intensity of the X-ray emission of galaxy groups scales with 
their mass, a large amount of them has not been detected at 
moderately high-z. Therefore we also expect a contribution 
to the overall signal from medium to low mass-groups at 
z>0. 2-0.3. 

Concerning the unresolved extragalactic CXB, it is difficult 
to distinguish its components with a simple spectral 
analysis mostly because of the poor energy resolution of 
X-ray sensitive CCDs (i.e. ~130 eV, for ACIS-I at 1.5 keV). 
However, cosmological sources leave a unique imprint in the 
power spectrum (PS) of the anisotropics of the fluctuations 
of the unresolved CXB in a way that is related to their 
clustering and volume emissivity properties. 
The unresolved CXB contains information about all those 
sources that have not been detected at the deepest fluxes 
reachable by deep surveys. Moreover, the amplitude of the 
PS is not only sensitive to the luminosity density of those 
sources, but it also provides information about their bias. 

In this paper we study the power spectrum of the unre- 
solved 0.5-2 keV CXB with the 4Ms CDFS data, the deep- 
est X-ray observation ever performed to date. We model the 
anisotropics of unresolved CXB with the state-of-the art re- 
sults on galaxy and AGN evolution models and observations 
as well as with modern hydrodynamical cosmological simu- 
lations. 

Throughout this paper we will adopt a concordance A-CDM 
cosmology with f2 m =0.3, Q.a=0-7, Ho=70 h km/s/Mpc and 
o"s=0.83. Unless otherwise stated, errors are quoted at la 
level and fluxes refer to the 0.5-2 keV band. We used as ref- 
er ence cumulative 0.5- 2 keV CXB flux the recent estimates 
of ILehmer et~aH (|2012T ), Scxs(0.5-2)=8.15±0.58xl0- 12 erg 
cm" 2 s" 1 deg" 2 . 



2 DATASET 

For the purposes of thi s work we used the 4 Ms Chandra 
exposure on the CDFS (|Xue et aUbOHl ). This is the deep- 
est Chandra observation ever performed and reaches a point 
source flux limit of ~10~ 17 erg cm -2 s" 1 . A detailed descrip- 
ti on of the data red uction and catalogue production is given 
in IXue et al] (|201lh . In order to improve the sensitivity on 
extended sources, we al so used the 3 Ms X MM-Newton ex- 
posure in the same field (|Ranalli et alj[2012l ). Since the scope 
of this paper is to measure the angular PS (see next section) 
of the purely diffuse, unresolved CXB, we masked the Chan- 
dra d ata from all the point sources detected by IXue et all 
(|201 ll ). The combination of the Chandra and XMM-Newton 
observations allowed us to remove, with unprecedented sen- 
sitivity, also extended sources down to 0.5-2 keV fluxes of 
2xl0~ 16 erg cm" 2 s" 1 (Finoguenov et al. , private com- 
munication). In order to employ the best-PSF area, we lim- 
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Figure 1. Left :The mask adopted for the ACF analysis in the observations of the CDFS. The white area is the resulting source-free 
area. Right : 0.5-2 keV image of the same area smoothed with a Gaussian filter with 3" width. 



ited our analysis to the inner 5' from the exposure weighted 
mean optical axis (a=03:32:27.316, £=-27:48:50.36). Point 
sources have been masked with circular regions of 5" ra- 
dius (which is large enough to mask >99% of the source 
flux at every off-axis angle investigated here) and extended 
sources have been masked out to R200 (i.e. radius within 
which the matter overdensity is >200). In this way the re- 
maining counts on the detector are made, with good approx- 
imation, by particle background and purely diffuse cosmic 
background only. With such a masking our image contains 
163940 counts, with an average occupation number of 1.7 
cts/pix. The total investigated masked area (80.3 arcmin 2 ) 
is shown in Fig. [T] together with the raw image of the same 
region. The source masking removes a fraction of the useful 
area of about ~25% and therefore the source-free area is of 
the order of 60 arcmin 2 . In principle the evaluation of the 
PS should not change if we increase the size of the mask. 
However by increasing the mask size we would run into an 
over-masking problem that would bias the estimate of the 
PS. If the masking were larger than ~30%-40% of the active 
pixels then the PS analysis w ould have been biased by the 
mask l|Kashlinskv et aL [2012). A smaller size of the mask 
would increase the overall power since it would be polluted 
by power from cluster outskirts and PSF tails of detected 
sources. In order to improve the signal to noise ratio of the 
weak signal we are looking for, we applied a careful filtering 
of the background. We have additionally filtered the events 
files provided by the CXC from particle flares by using the 
ciao routine Icjdean, which performs a 3<r clipping of the 
background light-curvfl Such a technique is more sensitive 
than the classical procedures adopted in standard pipelines. 
In order to reduce the noise introduced by low-count statis- 
tics, images were binned in pixels of 1.5" size. Images have 
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been produced in the 0.5-2 keV energy range, where Chan- 
dra has the largest collecting area and detector efficiency. 



3 EVALUATION OF THE POWER SPECTRUM 

In order to estimate the power spectru m we adopt 
the method d escrib ed in lArevalo et al.l ((2012) and 
IChurazov et al.l l|201lh that has been developed with the 
purpose of estimating the PS of 2D images with gaps. Here 
we give a brief description of the t echnique used to eva luate 
the PS, but we refer the reader to lArevalo et aL I (|2012h for 
a more comprehensive description. 

The first step is to derive the fluctuation field at any given 
angular scale 6, from which we derive a low-resolution PS 

by using a Mexican-Hat (MH) filter, F(x) = [1 - f^]e"^". 
The goal is to single out fluctuations at any given angu- 
lar scale and then compute the variance of the fluctuation 
field image, 8%{k). Our map has gaps introduced by fea- 
tures in the exposure maps and by the mask. The mask 
(M) is an image having value M=0 in the region of ex- 
cised sources and M=l elsewhere. The method corrects for 
gaps by representing the MH filter as a difference between 
two Gaussian filters with slightly different widths, convolv- 
ing the image and the mask with these filters and dividing 
the results b efore calculating the final filtered images. It has 
been shown (|Arevalo et al.|[2012h with numerical simulations 
that this method efficiently takes into account data masking 
also when a large fraction of the field is masked and, inde- 
pendently from the mask shape, the normalization and the 
shape of the PS is accurately estimated. At any given fre- 
quency k=l/d, we estimated the PS of the fluctuation field 
as follows: 

We define the fluctuation image at the scale 8 as: 
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where (G^o I, Gg 2 o I) represent the count-rate image (I) 
convolved with a gaussian filters with widths 0\^ (G) and 
M is the mas The width of the filter is chosen as 0i = 
Oy/1 + e and 02 = 0/V1 + e, where e «1. In this way the 
value in the brackets of eq[T]is dominated by fluctuations at 
the scale ~ 0. 

lArevalo et ail (|2012h demonstrated that PS of the fluctua- 
tion field at the frequency k—1/0 (Pa(k)) is related to the 
variance of image via: 



;a 2 (I e ) =< \5 s (k)\ 2 > . 



(2) 



e 2 7rfc 2 

In this work we used e =0.1, however a different value of 
e does not change the estimate of the PS, provided that 
e <<1. Note that in the calculation of the variance we 
have made no assumption on the distribution of the fluctu- 
ation amplitude, thus the method works also in the Poisson 
regime. This method is equivalent to the classical Fourier 
analysis; the only difference is that this treatment of the gaps 
and choice of the filters allow us to retrieve the underlying 
power-spectrum regardless of the shape of the mask. Unless 
otherwise stated, errors are given at the la level. In Fourier 
analysis, the uncertainty on power estimate is solely related 
to the accuracy with which one can estimate the variance. At 
every frequency, this depends on the number of independent 
elements in the Fourier space, and the intrinsic dispersion of 
the data. In our analysis, since we adopt a broad filter, we 
need to define, at every frequency k, the mean number of 
real Fourier elements used to evaluate the amplitude of the 
fluctuations < >. In principle the uncertainties can be 
approximated with op(k) ~ P2(fc)/v / 0-5 < Nt > ■ Following 
the extensive description of lArevalo et al.l (|2012h , within our 
formulation errors are computed with: 



a P (k) = 



y/f q (P(k)F? x dq N 



F 2 



(3) 



where Fkx is the Fourier transform of the filter, N is total 
number of pixels and N«=i is the number of unmasked pix- 
els. 

This filter-weigthed expected uncertainties is what we plot 
as error bars. 

Power spectrum points measured in this way are par- 
tially dependent on their neighbor frequencies as a result of 
the finite width of the power filter. Therefore, points can 
only be considered independent if they are measured at spa- 
tial scales that are sufficiently separated. In our case, where 
we used a factor of 2 separation in angular scale, it is appro- 
priate to use PS estimates and errors as independent points 
in the usual fitting methods. 



4 RESULTS 

4.1 Baseline method 

The power spectrum estimated with the technique described 
is Sect. [3] is shown in Fig. [2] We computed the probability 




The symbol o stands for convolution 



0(arcsec) 

Figure 2. Red squares : the raw PS of the source- masked 
fluctuations of the CXB as a function of the angular scale. 
Blue open circles : the (A-B)/2 PS. Black Triangles : The 
resulting Cosmic PS after subtracting the "spurious" power. 
Black symbols are artificially shifted for sake of clarity. 



that the observed data could be produced by a statistical 
fluctuation of a zero power signal by using a simple \ 2 test. 
We obtained a value of x 2 ~ 140 corresponding to a >10<7 
significance with respect to zero power. The high frequency 
part is dominated by a white noise component, while the 
signal increases towards lower frequencies (i.e. large scales). 
In Chandra images a large fraction of the fluctuation can be 
attributed to Poisson noise due to the low-count statistics. 

We estimat ed Pnoise directly from the data by using the 
A-B technique (Kashlins kv et aD [20051. For every pointing 
we sorted the events in order of time arrival and created two 
images, one including even events (A), and the other with 
odd events (B). We have then mosaiced all the A and B 
images. In this way we have two identical images with half 
the photons with respect to the full dataset, having both 
the same exposure, thus containing the same information 
on the signal but with their own noise. Any systematic er- 
rors should be manifested very similarly in the A and B 
images, and thus the (A-B)/2 difference images provide a 
useful means of characterizing the random (noise) proper- 
ties of the data sets. The PS of the difference image (A-B)/2 
is therefore nothing else than the PS of fluctuations due to 
Poisson noise. Source variability does not affect this proce- 
dure since the sampling of their light curve does not change 
significantly: the gap between two consecutive events is of 
the order of the read-out time of the detector, which is much 
shorter than the typical variability of X-ray sources. In Fig. 
[2] we show the PS of the Poisson noise compared with the 
PS of the raw image. As expected, the Poisson Noise shows 
a flat (white noise) spectrum. Throughout this paper we as- 
sume that both particle and galactic background contribute 
only to the Poisson noise since their flux is expected to be 
homogeneous at the scales sampled by this investigation. 
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The total cosmic PS (P 2 ,cxb) is therefore obtained with: 

P2,CXB = Pi,tot — Pi,{A-S)i%- (4) 

The estimated P2, cxb is plotted in Fig. [5] In the angular 
range sample here (3" < 9 < 1000"), the power increases 
with the scale, getting close to zero at smallest scale. The 
errors on P2, cxb have been obtained by propagating the 
errors of P 2 , to t and P 2 ,(a-b)/ 2 - 



5 MODELING THE CXB ANISOTROPIES 

In order to interpret the cosmic signal measured in the PS of 
the CXB, we used population synthesis models and recent 
observational results to predict amplitude and shape of its 
components. We want to test the assumption that the total 
observed fluctuations can be reproduced by a simple 3 source 
class model. Our hypothesis is that the amplitude and the 
shape of the PS are produced by: Shot Noise from individual 
undetected AGN and galaxies within the instrument beam 
(discrete source counts), clustering of galaxies and AGN be- 
low the limiting flux, clustering of undetected diffuse hot 
cosmological gas in large scale structures (IGM, i.e. unde- 
tected groups and WHIM filaments) . We then test the statis- 
tical robustness of our hypothesis with \ 2 test where every 
component is considered as a parameter. In addition we pro- 
vide an upper limit for the still undetected high-z sources. 
We assumed that, on the angular scales investigated here, 
fluctuations from the Milky Way diffuse emission does not 
contribute to the PS, i.e. this emission can be described as a 
constant flux on the detector and therefore contributing to 
the Poisson noise only. This hypothesis is supp orted by the 
PS m easurements of the soft CXB with ROSAT (|Sliwa et all 
l200lf ) that shows that the Galaxy produces signal only for 
the smallest harmonics of the PS. Moreover, cosmological 
sources (AGN, Galaxies and IGM) partially share the same 
environment, thus we include in the model a cross power 
term to model the amplification of fluctuations produced 
by their cross-correlation. In summary, the CXB PS can be 
modeled as: 



P2,CX B (k) = P 2 ,Siv(fe) + P 2 ,AGiv(fc) + 
P'2,GAL(k) + P 2 ,IGM(k) + P2,A,G,l(k), 



(5) 



where P2,SN(k), P 2 ,AGN(k), P 2 ,GAL(k) and P 2 ,iGAi(k) are 
the PS of shot noise and the clustering components of AGN, 
Galaxies and IGM, respectively. P 2 ,A,G,i(k) is the cross 
power term which contains the cross-correlation of AGN and 
galaxies, AGN and IGM, Galaxies and AGN, respectively. In 
the following sections we will describe the procedure adopted 
to model each spectral component. 

5.1 Preliminaries 

To decompose P 2l cxs in its primary componen ts we adopt 
the procedure described by iKashlinskvl (|2005h . Whenever 
a fraction of the sky of the order 9 <1 rad is concerned, 
one adopts the Cartesian formulation of the Fourier analy- 
sis. The fluctuation field of the CXB surface brightness (S) 
is defined as <5S=S(6 I ) — (5*). P 2 (fc) is related to the auto- 
correlation function C(9) through: 



where Jo(x) is the zero— th order cylindrical Bessel function. 
A very useful quantity is the mean square fluctuation within 
a finite beam of angular radius 9, or zero-lag correlation, 
which is related to the PS via: 

C(0) = {SS 2 ) e = 

1 f'°° (7) 
— P2(k)W TH (k9)kdk(xk 2 P2(k), 
2f Jo 

where Wth is the window function. Thus, in order to visu- 
alize the relative strength of the fluctuations at any given 
scale 1/k, a useful quantity is *Jk 2 P 2 (fc), that differs from 
the actual value of < 5s >e by a factor of the order unity 
that depends on the window function. 



5.2 Shot noise level and PSF modeling 

The shot noise is produced by discrete sources in the beam. 
The relative amplitude of this component scales as N -1 / 2 , 
where N is the mean number of sources entering the beam. 
P 2 ,Siv is therefore very important in surveys performed by 
instruments with high angular resolutions like Chandra. If 
sources are remov ed down to a flux Sn m , it can be shown 
l|Kashlinskvl 120051 ) that the shot noise component can be 
expressed as: 



2 dN x , c 
S 1^' 



(8) 



where d ^ is the differential logN-logS of all the X-ray point 
sources (i.e. AGN, Galaxies, Stars). In most cases, the sen- 
sitivity of the survey is not homogeneous across the field 
of view. This means that Su m in Eq. [8] is a function of 
the sky coordinates. To account for this selection effect we 
scaled the number counts by the selection function com- 
puted as follows: we computed the sensitivity map of the 
surveyed area b y imposing the sa me false source detection 
rate adopted bv lXue et al.l (|201 lh for catalogue production 
(i.e. P false <0.004). By using the masked maps described 
above we estimated the minimum count rate necessary not 
to exceed the threshold mentioned above. In Fig.[3]we show 
the selection function, in terms of fraction of field of view 
as a function of the flux limits. If r](S) is the selection func- 
tion for sources with flux S, the source counts of sources 
producing the shot noise is given by 



dN x 



a (1-77(5)) 



dN x 



(9) 



dS (Siv) dS 

a is a factor that accounts for the lower normalization ob- 
served in the CDFS with respect t o larger area fields an d 



with the iGilli et al l (|2007l ) mod el (ICappelluti et al.l | 2009h 
We measured a by co mparing thelLehmer et al.l ( 2012 ) 0.5-2 
keV logN-logS and the lGilli et all i|2007ft modeled logN-logS 
at the break flux and obtained a=0.72. We fixed the level 
of shot noise in our data by computing the integral in Eq. 
[H] by assuming that the only sources accounting for the SN 
are AGN and Galaxies. Therefore in this case we adopt: 



dN A 



dN x 



dS SN dS AGN 



+ 



dN x 



dS Gal 



(10) 



P 2 (k) = / C(9)J {k9)kdk, 



(6) 



dN/dS^Giv has been es t imate d from the population synthe- 
sis model of lGilli eT al, (2007) and dN/dSGai has been com- 
puted with the model of X-ray Galaxy evolution described 
in Sect. 15.41 In our calculations fluxes have been converted 



© 2002 RAS, MNRAS Q00.rTHT5l 



6 N. Cappelluti et al. 




0.001 1 1 — — 1 — — 1 — 

10-W 10-18 1Q -17 1Q-16 [0-15 

S (0.5-2) erg cm- 2 s"> 

Figure 3. Black Line The selection function for point sources in 
the inner 5' of the 4Ms CDFS. The grey shaded area represents 
the complementary of the selection function (l-rj(S)), highlighting 
the flux region from where the fluctuations are produced. 

into count-rates by assuming an average power-law spec- 
trum with spectral index T=1.7. With such an assumption 
count-rates are obtained by multiplying the flux by an en- 
ergy conversion factor ECF= 1.66xlO n cm 2 /erg. Note that 
a different choice of the spectral index of Ar ~ ±0.3 pro- 
duces a 5% variation of the ECF. The shot noise PS is scale 
independent (white noise) and since the signal from cosmo- 
logical sources increases with the scale, the contribution of 
P2,sjv to the total PS becomes more important at small an- 
gular scales. 

According to these prescriptions, we obtained an estimate of 
Psjv=1.58x 10~ 12 cts 2 s~ 2 deg -4 . The shot noise, as well as 
the other astronomical components of the PS, are affected 
by instrument PSF. This effect consists in a multiplicative 
factor that applies to the effective PS. In o rder to model it 
we use d the empirical approach proposed bv lChurazov et all 
l|201ll) where the PS of each modeled component is multi- 
plied at every frequency by factor: 

Ppsp{k) = [i + (fc/cWr' 

where k is the angular frequency. In Fig. [6] we show our es- 
timate of Psn, taking into account the effects of the PSF. 
When compared with the data we find an excellent agree- 
ment of our estimated Psn with the Cosmic PS at short 
scale (i.e. 8 <20") and we interpret this as a confirmation 
of our assumptions. Note that after subtracting the Psn 
from P2, tot (fc), the remaining signal can be completely at- 
tributed to clustered sources and therefore to the sources 
contributing to the unresolved extragalactic CXB. In order 
to evaluate the strength of our detection, we computed the 
probability that our data could be obtained with a random 
fluctuation of a null flat signal. If we subtract the SN compo- 
nent, at 9 >50", we still observe a strong signal. In order to 
determine the significance of the remaining signal, we fitted 
our data, on the whole angular range, with a P(fc)=0 model 



(zero power signal) to determine if our signal could be due to 
a statistical fluctuation of a null signal. With such a model 
we obtained x 2 /d-o.f.=43.5/13, thus rejecting our assump- 
tion at >4<r confidence level. Since such a signal increases 
with scale, we interpret this as clustered cosmic signal. 

5.3 Anisotropies from AGN clustering 

On the small angles sampled by our data, P2, cxb is re- 
lated to the unresolved CXB production rate AS/Az, and 
the evolving 3 D PS of the A GN, P3, AGiv(fc) via the Lim- 
ber's equation (|Peebleslll98"oh . At any given redshift the PS 
of AGN can be related to the PS of matter by knowing the ir 
redshift dependent linear biasing factor f b ( z ) . iKaiserl 1 1 9841 ) . 
Basically we have: 

^3,AGN (k,z) = b(z) 2 P 3M (k,z), (12) 

where Pa,M(k, z), is the matter 3D PS. P3,Af(fc , z) w as es- 
timated by using the CAMrfl l|Lewis fc Bridle! 120021 ) tool 
and including in the computation both the linear and th e 
non-linear components of the matter PS (|Smith et al.ll2003l ). 
Although a large amount of work on galaxies and AGN clus- 
tering is present in the literature, here we choose to combine 
this cosmological tool with bias measurements and predic- 
tions because of the poor sampling of the redshift space of 
observations. The 2D PS can be then obtained with: 

P2,AGN{k) = / ("4-) X 

JO \ dz J AGN ( 13 ) 

P 3 ,AGN(k[2nd A * (1 + gXTSf) dz 
c dt/dz [d A * (1 + z)] 2 1 + z' 

where d^ is the angular diameter distance, and the inte- 
gration is performed in the redshift range 0< z <7.5. A 
critical point of the analysis is the determination of the flux 
produced by undetected AGN at any given redshift. We de- 
rived ^ by usin g the CXB synthesis model published by 
iGilli et all (|2007l ). Briefly, the models take into account the 
observed luminosity function, fc- c orrect ions, absor ption dis- 
tribution and spectral shapes (see lGilli et alj|2007l . for more 
details) of AGN and returns the observed flux at any given 
redshift by using: 

JC r°° pz+dz rl it/ 

(14) 

where dt is the luminosity distance, L' is the luminosity 
measured in the 0.5(l+z)-2(l+z) keV range and dV/dz is 
the comoving volume element. The values of have been 
calcul ated with the interactive tool provided bv lGilli et al.l 
(2007j3. The population of AGN used to compute the 0.5- 
2 keV CXB production rate has the following properties: 
42< Log{L x ) <47 erg/s, 20<Log(A H ) <26 cm" 2 , and 
0< z <7.5. In addition, we included a high-redshift de- 
cline of the AGN s pace density (see e.g. iBrusa et al.ll2009l : 
ICivano et al.l l201ll ). They indeed modeled the evolution 
of AGN X-ray luminosity function (XLF) at z>2.7 with 
an exponential decay (0(L, z)=0(L,zo) x 10~ 0,43(z_zo) and 
zo=2.7) on top of the expected extrapolation from lower 

4 http://camb.info/ 

5 http:/ /www. bo. astro.it/~gilli/counts. html 
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Figure 4. CXB production rate for undetected AGN 
(Red continuos line). Our prediction for the undetected galaxies 
CXB production rate is shown as Green Line. The black line 
represents the emission from the IGM. 



redshift parametrization of lHasinger et all ([2005) . In Fig. [4] 
we show the undetected AGN CXB production rate as func- 
tion of the redshift. The bulk of their flux contributing to 
the unresolved CXB comes from z~l. 

Another important ingredient in the computation of 
P2, AGiv(fc) is the bias evolution of X-ray selected AGN. 
It has recently become clear that below z~3 AGN bias 
evolves like t hat of DM halos (DMH ) of mass of ~10 131 
Mq (see e.g. ICappelluti et al] |2012| . for a review) . How- 
ever, at higher redshift the bias factor of X-ray AGN is 
still unknown, while for optically selected QSOs this has 
been modeled up to z~5 (|Hopkins et al]|2007 ; iBonoli et al.l 
120091 ; Ida Angela et al]|2008l ) with qua dratic polynomials. I n 
Fig. [5]we compare the prediction of IBonoli et al.l {2009), 
rescaled to fit the z=0 X-ray selected AGN bias, and the 
bias computed fro m analytical models (jvan den Boschll2002l ; 
ISheth et al]|200ll ) for DMH with mass ~10 131 M . As one 
can notice, both curves fit in an excellent way the observa- 
tional data. Since the evolution at high-z of the AGN bias 
factor is unknown and still matter of debate, we assumed 
that the AGN bias evolves like the bias of DMH of mass 
~10 131 M Q up to z~10. The contribution of AGN to the 
overall PS is shown in Fig. [6] 



5.4 Anisotropies from X-ray galaxy clustering 

For anisotropies from galaxies we adopted a procedure sim- 
ilar to that used for AGN. The CXB production rate of 
galaxie^f] has been obtained by folding in Eq. [2]the z=0 lu- 



6 We considered "galaxies" all those X-ray sources with Ljf (0.5- 
2) kcV<10 42 erg/s. 



_ X-ray AGN (Cappelluti et al. 8012) 
15 _Ly-Break Galaxies (Ouohi et al. 2006) 
.Upper limits EBL (Robertson el al. 2010) 
-SFG (Adelberger el al. 2005) 



Figure 5. Black Continuous Line: Bias ev olution model 
adopt ed for AGN compared with data measured bv lAllevato et all 
l|201lh : (red dots). B l ack d ashed line: Bias evolution model for 
AGN of lBonoli et al] (I2009T) . Blue line: Bias evolution model as- 
sumed for X-ray emitting galaxies. Blac k pentagons : Bias mea - 
surements for star-forming galaxies from[Adclbcreer et al. (2005) 
and with blue ci rcles the me a surem ents of bias of high-z Ly- 
Brcak galaxies of lOuchi et al ] (|2004h . In green the upper-limit 
from EBL fluctuations of Robertson et al. (2010) 



minosity function measured bv lRanalli et al.l (|2005t ) with the 
evolution measured for star forming galaxies by Bouwens et 
al. (2011). The latter provides the most recent measurement 
of the evolution of star formation in the Universe up to z~10 
obtained with HST-WFC3. X-rays from non-active galaxies 
are mostly produced by low and high-mass X-ray binaries. 
Such objects with X-ray luminosities of > 10 37 erg/s can- 
not be detected individually in distant galaxies. The contri- 
butions from fainter discrete sources (including cataclysmic 
variables, active binaries, young stellar objects, and super- 
nova remnants) are well correlated with the star formation 
rate of the galaxy itself. However the ignition of X-ray activ- 
ity in the stellar population of galaxies has a delay from the 
burst in star formation of the order of the time scale of stel- 
lar evolution of the donor star in binary systems. Our model 
does not take into account effects of the delayed switch on 
of X-ray Binaries from the star formation as done by e.g 
Ptak et al. (2001); however, our representation allowed us 
to model the evolution of X-ray emitting galaxies up to z=10 
with the most recent results on st ar formation evolution. 
As shown by iRanalli et al.l (|2005l ) , the X-ray spectrum of 
galaxies can be represented by a simple power-law with 
spectral index V ~2. With such an approximation no k- 
correction is needed. 

In Fig. [4] we show the unresolved CXB production rate of 
galaxies as a function of the redshift. The galaxies contri- 
bution to the total flux of the unresolved CXB is dominant 
with respect to AGN at 0<z<6.5. Another parameter, that 
enters into the determination of the contribution of galax- 
ies to the X-ray PS, is their bias factor and its evolution. 
Most of the modern galaxy clustering analysis papers make 
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use of the Halo Occupation formalism and therefore it is 
not possible to derive an analytical formulation of the bias 
evolution. Moreover, galaxy clusterin g is complicated b y the 
luminosity dependence of clustering l|Ouchi et al.ll2004h . We 
therefore adopted a simple approximation for its determi- 
nation, by assuming that the comoving correlation length 
of star forming galaxies is consta nt, in the redshift range 
0<z<7.5, at the value measured bv lAdelberger et al.1 (2005) 
of ro= 4.5 Mpc/h and their spatial correlation function can 
be modeled with a power-law with 7=1.6. Within this sce- 
nario the bi as factor of ga laxies can be estimated at every 
redshift via l |Peebledll98Ch : 

b(z) = (Ts,G(z)/a8,r>M(z), (15) 

where <T8,g(z) is the rms fluctuations of the galaxies distri- 
bution over a sphere with radius of 8 Mpc/h and us,dm(z) is 
the same quantity for DM normalized to <78,.dm(z=0)=0.83. 
For such a power-law formalism it is possible to demonstrate 
that 

where J 2 =72/[(3- 7 )(4- 7 )(6-7)2 7 ] (|Peeblesl [l98ch ■ In Fig. [5] 
we compare our predicted bias evolutio n with the m e asure - 
ment of high-z Lyman break galaxies of lOuchi et al.l (|2004T ) 
and with upper limits derived by Robertson et al. (2010). 
As expected, galaxies are much less biased than AGN. 



5.5 Emission from Cosmological structures and 
WHIM 

In the inset of Fig. [5] we show that, after including the contri- 
bution of the components mentioned above, the power spec- 
trum still shows a prominent excess signal which, although 
with large error bars, increases up to a factor~3-4 toward 
low frequencies. Another significant contribution to the total 
signal of the CXB may arise from emission and clustering of 
unresolved galaxy clusters, groups and filaments. From the 
sensitivity maps of galaxy clusters we determined the flux 
limits for galaxy cluster/group detection and converted into 
a luminosity limit at every redshift. For galaxy clusters and 
groups luminosity , mass and temperature are related by 
sc aring relations. Th e relat ion adopted here are discussed 
in iFinoguenov et aTl l|2007T ). The luminosity limit can be 
then translated into a mass limit at every redshift. Thus 
our source detection ensures the removal of galaxy clusters 
and groups down to a Log(M)=12.5-13.5 M (i.e. kT<1.5 
keV). Thus only the low luminosity (low mass) and warm 
population of galaxy groups contributes to the unresolved 
CXB. Since this class of objects is very difficult to model 
analytically, we describe their properties using a set of mock 
maps from Roncarelli et al. (2012), who used a cosmologi- 
cal hydrodynamical simulation to define the expected X-ray 
surface brightness due to the large scale structures (LSS). 
The original hydrodynamical simulation (see the details in 
Tornatore et al. 2010) follows the evolution of a comoving 
volume of 37.5/i _1 Mpc 3 considering gravity, hydrodynam- 
ics, radiative cooling and a set of physical processes con- 
nected with the baryonic component, among which a chem- 
ical enrichment recipe that allows to follow the evolution of 
7 different metal species of the intergalactic medium (IGM). 




12 3 



z 

Figure 7. The Red line represents the lowest detectable mass in 
the 4Ms CDFS as a function of redshift obtained through scaling 
relations. 



From its outputs, Roncarelli et al. (2012) simulated 20 light- 
cones, covering the redshift interval < z < 1.5 with a size 
of ~0.25 deg 2 (roughly 4 times the size of the CDFs) each, 
and angular resolution of 3.5 arcsec. Each pixel of the maps 
contains information about the expected observed spectrum 
in the 0.3-2.0 keV band with an energy resolution of 50 eV. 
The emission coming from the IGM was computed assuming 
an emission from an optically thin collisionally-ionized gas 
(Apec in XSPEC) model and considering the abundances 
of the different metal species provided by the simulation. 
These maps/spectra have been convolved with the Chandra 
response in order to reproduce the effective Chandra count 
rates. Since our data are masked for galaxy clusters, we ap- 
plied to the simulations a source masking similar to that on 
the real data. The unresolved CXB production rate evalu- 
ated from simulations is shown in Fig. [3] We have simulated 
observations with the actual depth of the CDFS starting 
from the count rate maps described above, folded through 
a flat exposure map. We have added an artificial isotropic 
particle a nd cosmic background according t o the levels esti- 
mated bv lHickox fc Markevitcbl (|2006l .[2007). Random Pois- 
son noise was artificially added to the image and we ran 
a simple sliding cell detection with a signal to noise ratio 
threshold of 4. We have then excluded all the regions within 
which the overall encircled signal from sources is above 4 
sigma with respect to the background. Since galaxy clus- 
ters and groups are highly biased, in order to smooth out 
effects of sample variance we extracted the power spectrum 
from all the masked maps and averaged the results from 
all the realizations. Results of the PS modeling from unde- 
tected IGM is shown in Fig. [5] (magenta solid lines). A long 
standing debate in astrophysics is the possibility of detect- 
ing signal in emission or in absorption from the WHIM. In 
order to evaluate the contribution of WHIM to the overall 
PS signal of the unresolved CXB, we have extracted from 
our simulations the same lightcones but including only all 
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Figure 6. Left Panel. Red Triangles : Cosmic PS in the CDFS as a function of the scale. Red continuous line : PS of the shot noise 
(Galaxies and AGN). Green continuous line : PS of undetected AGN clustering. Cyan line: PS of undetected galaxies clustering. 
Magenta continuous line : Total power from undetected hot gas in the IGM. Blue — line : predicted clustering from of mini-quasars. 
The black line represent to overall model described in Eq. [5] The inset shows the data/model ratio as a function of the scale for a 
SN and AGN+Galaxy clustering model (Blue), SN+IGM clustering model (Red) and SN+ AGN, Galaxies and IGM clustering model 
(black). Right Panel: Amplitude of the fluctuations as a function of the angular scale with the same color coding as in the left panel 
with the addition of the contribution of IGM with 10 5 < T <10 7 K and overdensity S <1000 following the classical definition of WHIM 
Magenta dashed line to enhance its contribution to the total IGM flux. Yellow line : Cross-correlation contribution. 



those photons with a temperature 10 5 <kT<10 7 K coming 
from regions with matter overdensity S <1000. Such a se- 
lection is compliant to the classical definition of the WHIM, 
even if some denser clumps might be present also inside the 
filamentary structures (see the discussion in Roncarelli et al. 
2012). We have masked the WHIM photon maps with the 
same masks used for clusters, extracted the PS and averaged 
over the 20 realizations. The resulting modeled PS is shown 
in Fig. [6] (magenta dashed lines). 



5.6 Cross-Correlation terms 

As mentioned above, AGN, Galaxies and IGM partially 
share the same large scale structures. For this reason, fluctu- 
ations are boosted by their cross-correlation term. The 3-D 
cross-power spectrum (CPS) of the source populations 1 and 
2 is determined by: 

P 3l2 (k)=biblP 3 , M (k,z), (17) 

where bi and b2 are the bias factors relative to the source 
class 1 and 2, respectively. In analogy with the PS, the an- 
gular CPS of diffuse emission produced by two populations 
with background production rate ^ and j§ 2 , respectively, 
can be evaluated with the following form of the Limber's 
equation: 



P 3l2 {k[2-Kd A *{l + z)]-\z) dz 
c dt/dz [d A * (1 + z)] 2 1 + z' 

We have then computed P2 AGNJGM , P2 AGN Gal and 
P2 Gai IGM using the values of bias and emissivity described 
above. As far as IGM is concerned we adopted the bias evo- 
lution derived from simulations where biGM(z) = \/l + z. 
The overall cross-power term can be then expressed as: 

P2,A,G,I, P2 AGN IGM + P2AGN,Gal ^ 2 Gal,IGM • (19) 

The cross power term is plotted in yellow in Fig. [5] 
5.7 Accuracy of the PS model 

Our CXB PS model is shown in Fig.|S]where it is compared 
with data. Here we summarize the ingredients adopted to 
model the PS components and their limits. 

• At small separations most of the signal is due to shot 
noise produced by AGN and galaxies within the Chandra 
beam (red continuous line). 

• The AGN clustering is modeled by convolving the mat- 
ter PS for a concordance A-CDM cosmology, the AGN CXB 
production rate (dS/dz) derived from their X-ray luminos- 
ity function, and a linear bias evolution of DMH of mass 
Log(M)=13.1 Mq. Such a component increases with the 
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scale and its rms-fluctuations (right panel) peak on scales 
of the order of 40-50". Sources of uncertainties for such a 
component can be introduced by a) our assumption on the 
bias evolution that is unknown at z>2.5, b) the evolution 
of their luminosity function has been determined at high-z 
only for mos t luminous s ource s and c) the assumption at the 
basis of the iGilli et ail (|2007l ) population synthesis model. 
iFiore et al.l ( 20121 ') suggest that low- luminosity AGN may 
show a milder evolution at high-z with respect to observa- 
tions at high luminosity. This may result in an underestima- 
tion of the power at large scale. However, as discussed above, 
the peak of the undetected CXB flux is produced by AGN 
at z~1.5, where current observations are broadly consistent. 

• X-ray galaxies clustering (cyan line) is modeled by con- 
volving the matter PS for a concordance A-CDM cosmology, 
a model of the X-ray luminosity function obtained by as- 
suming an evolution similar to that of star forming galaxies 
starting from a measured z=0 XLF and a linear bias evo- 
lution of sources having a correlation length ro=4.5 Mpc/h. 
As for the AGN, such a component increases with scale and 
its rms-fluctuations peak (right panel) on scales of the or- 
der of 50". The sources of uncertainties introduced by our 
models are, like in the case of AGN, twofold and basically 
driven by our choice of the bias and XLF evolution. In the 
X-ray band the XLF of galaxies is known only up to z~1.5 
where it follows the evolution of star formation. The limit 
of our approach is on the adopted model describing their 
cosmological evolution, which is however based on the rea- 
sonable assumption that X-ray galaxies trace star formation 
in the Universe. The same argument applies to bias, where 
no measurements have been performed in the X-ray band. 
However, also in this case, observational proxies are consis- 
tent with our assumptions at the redshift where the bulk of 
the emission is produced, and therefore limiting the margin 



of uncertainty of our modeling. Moreover, as we describe be- 
low, our modeling of the shot noise of AGN and galaxies is in 
good agreement with the data, meaning that our estimated 
source counts below the flux limit is statistically robust. 

• The diffuse gas contribution (Magenta continuous line) 
is modeled by computing the PS of maps obtained by hydro- 
dynamical cosmological simulations which include a specific 
recipe of chemical enrichment. Sources detectable in real ob- 
servations have been excised from these simulations before 
computing the PS. Such a component is responsible for most 
of the large scale fluctuations that peak on a few arcmin- 
utes scale. For comparison, we also evaluated the expected 
contribution of WHIM (magenta dashed line) to the overall 
fluctuations and found that, on the largest angular scales, 
it contributes for up to ~l/3 of the total power. The main 
source of uncertainty for this component is definitely made 
by the actual metallicity of the IGM in the lowest density 
phase which has poor observational constrains. In addition, 
the largest scale clustering in the simulations may be un- 
derestimated because of the limited volume sampled in our 
cones. 

In order to test our statistical hypothesis, we performed 
a x 2 test and evaluated the improvement of the fit by 
adding one component after the other. We considered ev- 
ery component (which have their amplitude and shape fixed 
by our assumptions) as if they were a parameter of the 
fit for the evaluation of the degrees of freedom. The first 
test has been performed to evaluate if the observed fluctua- 
tions could be explained by the shot noise component only. 
As discussed above, such an hypothesis is rejected at >4<r 
confidence level. Then we tested if, together with a shot- 
noise component, clustering of AGN and galaxies (including 
their cross power) improved the fit. In this way we obtained 
X 2 /d.o.f.=17.5/12. We also performed an F-test to evaluate 
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Model 


X 2 /d.o.f 


. P(F-test) 




(1) 


(2) 


PSN 


43.4/13 




P,s n+PgcI+Pagn 


17.5/12 


1.2X10 -3 


PsN+PlGM 


18.3/12 


1.6xl0" 3 


Psn+PgcI+Pagn+Pigm 


7.8/11 


2.7xl0" 3 * 


PsN + PGal+PAGN+PlGM+Pmq 


8.3/10 


* 



Table 1. (1) The x 2 /d.o.f. value for model with different components include and (2), the null hypothesis probabilty obtained by adding 
an additional component with respect to the model with n-1 d.o.f. *:P(F-test)=7.9x 10 — 5 with respect to the SN only model. 



the probability that the obtained x 2 could be obtained by 
a statistical fluctuation of the SN only model and obtained 
P(F-test)=1.2x 10 -3 , thus the inclusion of such a component 
significantly improves the fit. However, the relatively high 
X 2 /d.o.f. value suggests that additional or different compo- 
nents are required. 

A visual inspection suggests that the strongest feature in 
the PS is produced by the IGM feature. Thus we fitted the 
data with SN+IGM model and obtained x 2 /d.o.f.=18.3/12 
corresponding to P(F-test) = 1.6x 10 -3 with respect to the 
SN only model. For such a model we show the data/model 
ratio in red in the inset of Fig. [6] 

By adding the point source clustering (and all the cross 
power terms) to the latter model we obtain x 2 /d.o.f.=7.8/ll 
and therefore the F-test, computed to the SN+IGM model 
provides P(F-test)=2.7x 10 -3 , thus providing a further ~3<r 
improvement of the fit. If we compare this fit with the SN 
only model the F-test probility is P(F-test)=7.9x 1CT 5 , thus 
providing a significant improvement of the overall fit quality. 
The overall data/model ratio is shown in black in the inset 
of Fig. H 

According to the values reported in the table, we can safely 
confirm that the unresolved CXB in the CDFS can be ex- 
plained with random and clustered signal from undetected 
point sources (AGN and galaxies), IGM clustering compo- 
nents plus a cross correlation term. 



6 FLUCTUATIONS FROM EARLY BLACK 
HOLES? 

Potential contribution to the unresolved CXB and its struc- 
ture can also come from very high redshifts overlapping with 
epochs usually identified with first stars era. It is widely ex- 
pected that fragmentation within the first collapsing proto- 
galaxies was much less efficient so that first stars were sig- 
nificantly more massive and short-lived and could have left 
behind non- negligible abundance of accreting black holes. 
Since, in addition, the first black holes could have also 
formed directly during the first stars era, these populations 
may supply an additional, and potentially measurable, con- 
tribution to the unresolved CXB and its fluctuations. 

Spitzer-based studies have revealed significant levels of 
source-subtracted cosmic infrared background (CIB) fiuctu- 
atio ns which were proposed to originate in the first stars 
era IjKashlinskv et alj 120051 '). Indeed the remaining known 
galaxy populations have been shown t o produce significantly 
lower levels of the CIB fluctuations (jHelgason et al.l 120121 ) 
and there is no correlation in the large-scale structures be- 



tween the Spitzer CIB maps at 3.6 and 4.5 micron and 
HST/ACS likely pointing to the high-z origin of the ex- 
cess CIB |Kashlinskv et all 120071 1 . The excess fluctuations 
have been confirme d in the AKARI based analysis extend- 
ing to 2.4 micron {Matsumo to et all l201lT ) and the signal 
has now been measured to extend to ~ 1° exceeding the 
power over the rem aining normal ga laxies by well over an 
order of magnitude (|Kashlinskv||2005l ) . The level of the fluc- 
tuations at 3.6 micron is SFcib ~ 0.07 nW/m 2 /sr and their 
energy spectrum appears to follow the Rayleigh- Jeans law, 
SFcxb oc A -3 , between 2.4 and 4.5 micron. If the excess 
CIB fluctuation arises at high z, the sources producing it 
would have projected angular density of ~ (1 — 3)arcsec -2 
(Kashlinsky et al 2007). Under the assumption that the 
high-z sources measured bv lKashlinskv et alj (|2012l ) rapidly 
evolve into miniquasars ( or early black holes) , we estimated 
that, if the CIB/CXB flux ratio is constant along the cos- 
mic time, with 3.6 extr agalactic CIB flux of ~6 nW 
m~ 2 sr" 1 (|Kashlinskv| 120051 ) and ex tragalactic CXB flu x 
of 8.15xl0" 12 erg cm -2 s" 1 deg" 2 i|Lehmer et all 120121 ) . 
we obtain |£J_b ^ 2 20. Thus we evaluated the contribu- 
tion to the C&B of the CIB fluctuations to be of the order 
~3xl0~ 13 /(i+z) erg cm" 2 s' 1 deg" 2 (i.e. ~0.5% of the 
CXB if these sources have z~7.5-15fl. We have therefore 
predicted the expected PS of these sources by folding their 
CXB production rate in Eq. [13] and integrating in the red- 
shift range 7.5-15. 

The result of this prediction is shown with a blue line 
in Fig. El A x 2 fit of our 5 components model provides 
X 2 /d.o.f=8.3/10. Although the inclusion of such a compo- 
nent does not improve the quality of the fit, the statistics 
does not allow us to reject the hypothesis that these sources 
contribute to the observed fluctuations. These sources would 
contribute very weakly to the shot noise, and the upper limit 
of their contribution is not larger than the uncertainty on 
the measured PS. Thus at la we have PsN,mq <1 X 10~ 13 cts 2 
s -2 deg _4 (i.e. the mean value of the uncertainty on those 
scales). According to EqUl we can estimate the source den- 
sity of putative miniquasars with Psjv /^hm-n ~N(> S m i n ). If 
Smm = 10 _20 erg cm -2 s _1 , where most CXB models predict 
the saturation of the source counts, we estimate N<60.000 
deg~ 2 . If these sources shine in the redshift range 7.5-15 
ijKashlinskv et al.ll2012h . then their comoving volume den- 
sity is ^9 xl0~ 5 Mpc~ 3 . For comparison in the same redshift 
range the iGilli et alj (|2007l ) model predicts, for X-ray se- 
lected AGN, in the luminosity range sampled by our fluctu- 
ations (Log(Lx)<43.58 erg/s, at z=7.5-15, for sources with 

7 Assuming b(z)=\/l + z 
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observed -20^LogS(0. 5-2)^-17 erg cm 2 s 1 ), source den- 
sities of the order 1.5xl0" 4 Mpc" 3 and 8.4xl0~ 5 Mpc -3 
in the case of flat and declining evolution, respectively. The 
upper limit derived above suggests that the declining evo- 
lution of AGN is a good representation of the evolutionary 
track of these remote sources. 



a CXB flux produced by undetected AGN of ~2.0xl0~ 13 
erg cm -2 s _1 deg -2 . 

We also tried to probe different evolution scenarios but 
our data do not allow to significantly constrain the behavior 
of the AGN XLF at high-z. In Fig. |9]we show the predicted 
LogN-LogS that, according to our model, satisfies the ob- 
served fluctuations compared with recent observations. 



7 CONTRIBUTION OF UNDETECTED 
SOURCE POPULATIONS TO THE CXB 

By interpreting the observed behavior of the unresolved 
CXB fluctuations, we have developed a model which is able 
to explain the nature of the unresolved CXB. The fluctu- 
ations observed here are reproduced in the PS if the clus- 
tering recipes described in the text are combined with the 
unresolved CXB production rates shown in Table [2] Our re- 
sults show that, in the 0.5-2 keV band, the effective fraction 
of the unresolved extragalactic background is of the order of 
12% of the total. In Table [2] we show the CXB flux that our 
model predicts to explain the fluctuations together with the 
fractions of total and unresolved CXB produced by every 
undetected source population. As one can notice in Table 
[2] the bulk (~56%) of the unresolved CXB flux is made by 
unresolved clusters, groups and the WHIM which accounts, 
by itself, for ~17% of the unresolved flux. 
For point sources our model predicts that AGN and galaxies 
contribute, all together, for the remaining flux of the unre- 
solved CXB, with galaxies and AGN producing ~25% and 
~20% of the unresolved flux, respectively. 
Moreover, our data cannot exclude that a sizeable fraction 
of the unresolved CXB could be produced by a population 
of still undetected high-z sources, likely black hole seeds. 



8 DISCUSSION 

In this paper we presented the measurement of the angular 
PS of the fluctuations of the unresolved CXB in the 4Ms ob- 
servation of the CDFS in the angular range < 10'. Poisson 
noise and spurious signals have been modeled and removed 
from the measured PS. We performed a spectral decompo- 
sition analysis and showed that after removing the low fre- 
quency signal, which can be attributed to the shot noise of 
unresolved sources that randomly enter the beam, the am- 
plitude of the fluctuations with extragalactic origin account 
for ~12.3% of the CXB and the significance of the detection 
of these cosmic fluctuations is >10cr. In the next section we 
briefly discuss to properties of the populations producing 
the unresolved CXB fluctuations. 



8.1 The population of undetected AGN 

For AGN we folded the observed evolution s c enari os with 
the population synthesis model of iGilli et al.l l|2007T l and a 
simple recipe for bias where AGN are tracing DMH with 
mass Log(M) = 13.1 M . 

This population of AGN has a space density that ex- 
ponentially declines above z=2.7. The CXB production rate 
necessary to produce the modeled PS yields to a fraction of 
the unresolved 0.5-2 keV CXB flux of ~19%. We predicted 



8.2 The population of undetected X-ray galaxies 

Galaxies are the most numerous population of objects con- 
tributing to the unresolved CXB (see Fig. the power 
produced by such a population is lower that that of AGN 
since they are less biased even if they produce more CXB 
flux. In the soft X-rays, galaxies have been observed up to 
z~l (see e.g. lLehmer et al]|2007l ). and therefore their high-z 
space density is unknown. We have developed a toy model 
for the galaxies XLF where their evolution follows that of 
the star formation in the universe. With this model we es- 
timated that galaxies contribu te to ^25% of t he un resolved 
CXB flux. In a recent paper, iDiikstra etail (|2012h used a 
modeling similar to ours and found that, in principle, X-ray 
galaxies could produce all of the unresolved 1-2 keV CXB. 
However they did not consider the contribution of other un- 
detected sources that we have shown to produce a large frac- 
tion of the unresolved CXB. Overall, the predicted source 
counts of AGN a nd galaxies are in goo d agreement w i th the 
measurements of lXue et alj (|201 lh and lLehmer et all (|2012T ) 
in the same field. 



8.3 CXB from IGM and WHIM 

On scales larger than 100" our analysis shows that the main 
contribution to the unresolved CXB fluctuations is due to 
IGM emission produced by undetected groups, clusters and 
the WHIM. This has been estimated with cosmological hy- 
drodynamical simulations that, together with structure for- 
mation, include feedback mechanisms that pollute the IGM 
with metals which are responsible for the X-ray emission 
l|Roncarelli et al.ll2012T ). Such a component produces ~50% 
of the unresolved CXB. 

We determined that ~ 1/3 of the IGM flux is produced by the 
WHIM, where the so called "missing baryons" are expected 
to lie. However the WHIM definition taken from simulations 
already applies a density cut (S < 1000) that is meant to 
mimic roughly the effect of excising detected sources. In re- 
ality, dense clumps may be present outside virialized objects 
so we must consider our determination as a lower limit of the 
total WHIM contribution (see the discussions in Roncarelli 
et al. 2006, 2012). 

Our model thus predicts that the flux of the WHIM in 
the CDFS is of the order 1.7xl0~ 13 erg cm" 2 s" 1 deg" 2 (i.e. 
2.3% of the total CXB flux) and produces a signal peaking 
on a few arcmin scale. Such an estimate is a bout one order 
of ma gnitude lower than what measured bv lGaleazzi et all 
(2009 ) (i.e. 12% of the overall diffuse emission in shallower 
0.2-0.4 keV XMM- Newton observations) where they did not 
model the contamination from undetected sources. However 
such an estimate relies on the output of the simulation and 
is very sensitive to t he me tallicity of the WHIM. According 
to iRoncarelli et all (|2012T l, different recipes for the metal 
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Table 2. (1) CXB production rate of every class of undetected sources. (2) The percentage of the overall and (3) unresolved CXB 
produced by every population. In the last line we show the cumulative values, the ranges refer to Model I and II. 



Component 


S 

xl0 -13 erg cm" 2 s -1 deg -2 
(1) 


% CXB a 
% 
(2) 


% Unres. CXB 

% 
(3) 


AGN 


1.97 


2.4±0.2 


19.3±1.3 


Galaxies 


2.51 


3.1±0.2 


24.6±1.7 


IGM 


5.70 


7.0±0.5 


55.9±3.9 


''WHIM 


1.70 


2.1±0.1 


16.7±1.2 


Early-BH 


<0.35 


<0.5 


<3.4 


Total 


10.18 


12.4±0.9 


100 



Computed using a total 0.5-2 keV CXB flux of 8.15±0.58x 10 -12 erg cm" 2 s" 

b The WHIM flux is included in the IGM flux. 



deg -2 jLehmer et al.ll2012Tl . 




10-" 10-' 6 
S(0.5-2 keV) (erg cm- 



Figure 9. The 0.5-2 keV logN-logS used to reproduce the ob- 
served unresolved CXB fluctuations. In green AGN, blue galax- 
ies and, in black — shaded, the sum of the AGN and galax- 
ies, the uncertainty is attributed to the count-rate to flux con- 
version. The red s haded area is the logN-logS measured by 
iLehmer et"aD \20l± in the CDFS. The red line represents the up- 
per limit of the logN-logS of miniquasars. The cyan shaded area 
represents the expected counts from the fluctuation analysis of 
iMivaii fc Griffiths! i2002T) . In the inset we show a zoom onto the 
miniquasar upper-limit logN-logS. The blue triangle is the ex- 
pected number of X-ray sources, at the flux limit of the 4Ms 
CDFS, produced at z<10 from collapse of POPIII stars accret- 
ing at the Eddington Limit. The red triangle is number of X-ray 
sources in case of of direct collapse (Quasi-stars) accreting at 



enrichment of the WHIM may lead to a variation up to a 
factor 3 in overall emissivity of the WHIM. More information 
on such a component of the Universe will be possible if, 
for example, contamination of undetected cluster could be 
excised from X-ray maps by masking also optically detected 
groups with mass Log(M)>12-12.5 Mq. 



8.4 Very high-z sources 

Finally, we speculated on the possible existence of a 
population of high-z miniquasars (or early black holes) 
born from the collapse of early massive objects. Although 
we did not detect their signature, we placed an upper limit 
to their contribution to the CXB (<3x 10 -13 /(l+z) erg 
cm -2 s -1 deg -2 , z>7.5). We estimated that these sources 
would follow the declining evolutionary track of AGN with 
Log(L x )<43.58 erg/s. 

Our observations are not sensitive to the faint fluxes 
expected from these sources, and thus we could only place 
upper limits. However, since such a population peaks at 
z> 7.5, their fluctuations should peak on scales of the order 
of several tens of arcminutes, where the contribution from 
shot-noise is relatively weak also in shallower surveys. On 
such a scale, foreground source population PS significantly 
dims and therefore their detection would be possible. To 
conclude, we determined that at fluxes of the order 10 -20 
erg cm -2 s -1 , their number density is order of ~60000 
deg -2 which means that, under the assumption of Euclidean 
logN-logS, at the flux limit of the CDFS their number 

density is of the order 1-2 deg -2 . 

Using the theoreti cal predictions of IVolonterU (|2010T ). 
iTreister et alj (|201ll ) computed the expected number den- 
sity of the first X-ray sources at z~7-10 at the 4Ms CDFS 
flux limit. In the case of Population III stars remnants 
accreting at the Eddington limit they find 4.5 deg~ 2 , while 
in th e case of direct collapse (quasi-stars, IBegelman et al.l 
2008) with accretion regime at A^dd^O.3 their expected 
number density is 1.9 deg -2 . In Fig.|9]we show a comparison 
of these predictions with our upper-limit which, at the 
zero-th order, is favoring the direct collapse scenario. 

Our work suggests that future deeper observations on 
wider fields would allow us to improve the sensitivity of the 
PS measurement by reducing the Shot-Noise and masking 
fainter clusters. Moreover a deep survey with a flux limit 
comparable or deeper to that of the 4Ms CDFS covering 
1-2 deg 2 , would allow us a direct detection of the WHIM 
feature and, by measuring the large scale shape of the PS, 
investigate the very high-z X-ray Universe. A possible future 
X-ray m ission like the propos ed Wide Field X-ray Telescope 
fWFXT [Giacconi et al.ll2009T ) will be able to study the PS 
of the unresolved CXB with high precision given the large 
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collecting area (low photon noise) and the large field of view 
(low cosmic variance). 
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